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ABSTRACT 

Context. With growing evidence for the existence of very massive stars at subsolar metallicity, there is an increased need for corre¬ 
sponding stellar evolution models. 

Aims. We present a dense model grid with a tailored input chemical composition appropriate for the Large Magellanic Cloud. 
Methods. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momen¬ 
tum by magnetic fields, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial 
masses from 70 to 500 Mg and with initial surface rotational velocities from 0 to 550 km/s, covering the core-hydrogen burning phase 
of evolution. 

Results. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homo¬ 
geneous evolution occurring for the fastest rotating models. Above 160Mo, homogeneous evolution is also established through mass 
loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment 
is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above ~ 150Mo at 
zero age, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proxim¬ 
ity of these models to their Eddington limit. This leads to a maximum zero-age main sequence surface temperature of ~ 56 000 K, 
at ~ 180 Afo, and to an evolution of stars in the mass range 50Mo... lOOAfo to the regime of luminous blue variables in the HR 
[spell out] diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically 
homogeneous evolution of stars above ~ 180 Mq. 

Conclusions. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents 
stars at LMC metallicity from evolving into pair-instability supemovae. The corresponding spin-down will also prevent very massive 
LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses. 
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1. Introduction 


Massive stars, wit h initial masses abo ve ~ SM©, are power¬ 
ful cosmic engines (iBresolin et al.ll200^ . They produce copious 
amounts of ionizing photons, strong stellar winds, energetic fi¬ 
nal explosions, and most of the heavy elements in the Universe. 
The most massive amongst them conduct a life very close to 
their Eddington limit, and may thus be prone to become unsta¬ 
ble. They are thought to be able to produce the most spectacular 
stellar explosions, like hypernovae, pair-instability supernovae, 
and long-duration gamma-ray bursts (lLangeill2012h . 

While the value of the upper mass limit of stars is presently 
uncertain (ISchneider et al.ll2014ll . there is a great deal of evi¬ 
dence of stars with initial masses well above lOOM© in the local 
Universe. A number of close binary stars h ave been found with 
component initial mass e s above lOOM p ( Schnurr et al.l l2008l 
l200^[Sana et al.l|2013bh . lCrowther et al.ld^lOl) proposed initial 


masses of up to SOOMq for several stars in the Large Magellanic 
Cloud (LMC) based on their luminosities. 

We present stellar evolution models of very massive ro¬ 
tating core-hydrogen burning stars, with initial masses up to 
500 Mq. These models will be used for compar ison with the VET 
ELAMES Tarantula survey (lEvans et al.ll201 Ih . which observed 
more than 800 O and early-B stars in the 30 Doradus region lo¬ 
cated in the Large Magellanic Cloud, to study the effects of ro¬ 
tational mixing and mass loss on the evolution of very massive 
stars. 

Stellar evolution models f or very massive stars have already 
been presented (see Table 1 in iMaeder & Mevn^(l201 ih and ref¬ 
erences th erein). In particular, we refer to the rotat ing stellar 
models of ICrowther et al.l (l2010l) : lYusof et ^ (1201 3h . who cal¬ 
culated models with initial masses up to 500 Mq. It is the aim of 
this paper to present a grid of stellar evolution models that has a 
dense spacing in mass and rotation rates, such that it is suitable 
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for forthcoming population synthesis calculations. In this sense, 
our mo dels are an extension of the LMC models of iBrott et al.l 
(1201 lah to higher masses. In the mass range considered here, 
stellar wind mass loss and the proximity to the Eddington limit 
play prominent roles. 

In Sect. |2]we briefly describe the employed stellar evolution 
code, including the input parameters for our calculations. Fol¬ 
lowing that, we present and discuss the stellar evolution mod¬ 
els in Sect. [S] including the evolution of their mass and their 
surface abundances, with emphasis on chemically homogeneous 
evolution, and we compare our models with previous work in 
Sect. |4] Section |5] contains our short summary. In Appendix lAl 
we present isochrones derived from our models, as functions of 
age and initial stellar rotation rate. 


2. Input physics and assumptions 


iBrott et ^ (1201 lal) computed three grids of stellar evolution 
models for different metallicities (Galactic, SMC, LMC) for ini¬ 
tial masses up to 60 M©, which were compare d with the results of 
the VLT-FLAMES Survey of Massive Stars (lEvans et al.ll2005l 
l2006h . In particular, the convective core overshooting parame¬ 
ter and the e fficiency pa r ameter s for rotationally induced mix¬ 
ing used by IBrott et ^ (1201 lal) were calibr ated to reproduce 
the results of this survey ( Hunter et al.ll200^ . The dense spac- 
ing in initial mass and rotational velocity used for the grid of 
IBrott et ^ (1201 lal) opened the door for sta t istical t ests of the the¬ 
ory of massive star evolution (IBrott et al.l (1201 Ibl) : S chneider et 
al., in preparation). Here, we extend the LMC grid of iBrott et al.l 
up to initial masses of 500 Mq, using the same numerical code, 
physical assumptions, and a similarly dense grid spacing. 

To calculate rotating stellar evolution models, we use our 
one-dimensional hydrodyna mic binary stellar e yolution code 
(BEC) whi ch is described in iHeger et al.l (l2000l) : iPetroyic et al.l 
((2005 ), and lYoon et al.l(l2012l) . It contains a detailed treatment of 
rotation, angular momentum transport due to internal magnetic 
fields, and stellar wind mass loss. Tfie code solyes all fiye stellar 
structure equations throughout the stellar interior, including the 
stellar enyelope up to a Rosseland optical depth of r = 2/3. Con- 
yection is considered throughout the star using the non-adiabatic 
mixing-length theory. Our code is suited to treating stars close to 
the Eddington limit, and to descri be the effects of enyelope infla¬ 
tion w hich occur in this situation (llshii et alJl999l : IPetroyic et al.l 
I2nn6h . 


2.1. Chemical composition 

The initial chemical composition for our models is chosen ac¬ 
cording to corresponding obserya tions of young ma ssiye stars 
and of Hii-regions in the LMC (IBrott et al.l l201 lah . We thus 
adopt initial mass fractions for hydrogen, helium, and the sum 
of all metals of Y = 0.7391, Y = 0.2562, and Z = 0.0047, 
respectiyely, wi th a non-so lar metal abundance pattern (see Ta¬ 
bles 1 and 2 in IBrott et al.l) . We note that the applied opacities 
and mass loss rates (see below) are scaled with the LMC iron 
abundance, not with the total metallicity, which is reduced by 
0.35 dex with respect to that of the Sun. 

2.2. Convection and rotational mixing 

Conyection with a m i xing-length pa rameter of aMur = 1-5 
(iBohm-Vitensel 119581: lLangeiiTl991l) is applied as well as 
semi-conyection with an efficiency parameter of crsEM = 1 


(iLangeret alJll98^ iLangeil [19911) . In addition to conyection, 
conyectiye core oyershooting is included wi th anvRr = 0 . 335 lo- 
cal pressure scale heights, as calibrated in iBrotte^^J 20]Ta) 
with the rotation al properties of B-type stars ( Hunter et al.ll2008 ; 

I Vink et aLll2010l) from the VLT-FLAMES suryey. While no ob- 
seryational calibration of the oyershooting parameter exists for 
stars of the considered mass range, we point out that the role of 
oyershooting in our models is minor because of the large con¬ 
yectiye core mass fract ions of yery massi ye stars. 

Rotational mixing (iHeger et al.ll2000l) is consider ed with the 
efficien cy parameters of fc — 0.0228 and ^ = 0.1 (IBrott et al.l 
1201 lal) . The most significant process causing rotationally in¬ 
duced mixing in our models is the Eddington Sweet circulation. 
Furthermore, the transport of angular momentum b y magnetic 
fields due to the Spruit-Taylor dynamo (ISnrui I 2 OO 2 I) is applied, 
which is assumed h ere not to lea d to an additional transport of 
chemical elements (ISpruid l2006h . Since the magnetic torques 
lead to a shallow angular yelocity profile in our models, the ef¬ 
fects of the shear instability, although included, are quite limited 
during the main sequence eyolution. 


2.3. Mass loss 

The eyolution of yery massiye stars is intimately connected to 
their mass-loss behayiour. Mass loss of yery massiye stars, for 
which few obseryational constraints exist, is a yery actiye field 
of research. Here, we proyide a brief oyeryiew of the mass loss 
prescriptions adopted in our calculations. 

From the zero-age main sequence up to a surf ace helium 
fraction Yg - 0.4, we use the mass-loss predictions by lVink et al.l 
(1200(1 12 OOII) for O and B stars. These rates are yalid for stars 
of one million solar luminosities or less, i.e. stars below 80 Mq 
( iMokiem et al.ll2()(y^ . which are not yery close to their Edding¬ 
ton limit (i.e. haye E £ 0.3; see Eq. El). Empirical tests of 
these predictions depend critically on the presence and proper- 
ties of small scale str uctures in the outflows, known as clumping. 
iMokiem et al.l (l2007l) showed on the basis of Her and He ii 44686 
analysis that the predicted rates agree with obseryations if the 
material is concentrated in clumps that haye a h-A times higher 
density than in a smooth outflow. Other analyses, which include 
modelling of ultrayiolet resonance l ines, de r iye clumping factors 
that may reach yalues of 10 (e.g. i Boured (l2004l) : iBouret et al.l 
(l2005ll2()T^ : lFullerton et akl (l2006l) l. In this last case, the Vink 
et al. prescription adopted here may oyerestimate M by about a 
factor of 2. An improyed hydrodynamical treatment shows that 
for normal O stars the mass loss rates ma y be somewhat lower 
(iMiiller & Vinkll200llMuiires et al.ll20T2l) . 

Predic tions for stars in t he 40-300 Mq range haye been pre¬ 
sented by IVink et akl (1201 ll) . Objects at the upper mass end of 
this range may haye a yery high Eddington factor E. It has been 
found that, at solar metallicity, the wind strengths agree with the 
standard Vink et al. recipe for objects that haye an Eddington 
factor for Thomson scattering of Eg ^ 0.7. For higher Eg yalues 
these new predictions show an upturn in the mass-loss rate, lead¬ 
ing to rates that are higher by up to about a factor of 2 compared 
to the Vink et al. yalues used here. The authors associate the up¬ 
turn with the stellar wind becoming optically thick, i.e. leading 
to spectral morphology which is typical for Wolf-Rayet stars of 
nitrogen subclass (WN) showing hydrogen emission lines. Inter¬ 
estingly, t his predicted upturn may haye been confirmed obser- 
yationally dBestenlehner et a02014l) . The relation between mass 
loss and Eddington factor has also been explored for l a te-WN 
stars by iGrafener & HamannI (l2008h and iGrkfener et alJ ( 2014^ 
They find a behayiour that is similar to the results of Vink et aP 
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(1201 ih . though t he onset of Wolf-Ray et type outflows occurs at 
lower Fe values. iGrafener & HamannI however report a tempera- 
tu re dependence of their mass-loss rates that is steeper than that 
of IVinket al.](l201lh . 

Since the mass loss predictions for large Eddington factors 
cannot yet be implemented unambi guously in t o stell ar evolu¬ 
tion calculations, we extrapolat e the ^nket^j200]J) rates for 
stars above IO^Lq. We note that ICrowther et alT(l201 Oh found the 
Vink et al. rates to agree within error bars with those observed 
in stars of up to ~ 10^L© found in 30 Doradus. For objects 
in the range 60... 100 M© the objects are close to the model- 
independent mass-loss transition point between optically thin 
and thick winds. In this range mass-loss rates h ave recently been 
calibr ated with an uncertainty of only ~30% (IVink & Grafened 

IMl- 


The Vink et al. mass-loss prescription shows a bi-stability 
jump at about 25 000 K, leading to an increase of the mass-loss 
rate by a factor of 5 for stars of spectral type B1.5 or later. We 
include this bi-stability jum p in our calculations (cf. Brott et aP 
(1201 lah l. Additionally the iNieuwenhuiizen & de Jagej (Il990l) 
empirical mass-loss rate is applied to cope with an in¬ 
crease in mass loss when approaching the H umphreys- 
David s on lim it (HD limit). The transition from tVink et alJ 
(l2000i 1200 ih to INieuwenhuiizen & de .lageil (Il990h occurs at 
any effective temperature s maller than 22000 K where the 
Nieuwenhuiizen & de Jagej (1 19901) mass-loss rate exceeds the 


Vink et al J dioooi 2001 ) mass-loss rate. 


To account for Wolf-Rayet mass loss, the iHamann et al.l 
(1 19951) mass-loss rate divided by a factor of 10 is applied for 
the su rface helium fraction > 0.7. Figure 1 in lYoon et ^ 
(1201 Ol) shows that this corresponds well to the W olf-Rayet mass- 
loss rate proposed by iNugis & LamersI (l2000l) for Wolf-Rayet 
masses in the range 5Mq to 2QMq. For surface helium mass 
fractions from 0.4 < Tj < 0. 7 the mass - loss r a te is l inearly 
interpolated between either the Vink et aP (l2000l 1200 ll) or the 
Nieuwen huiizen & de .lageil (Il990h mass-loss rate and that of 


Hamann et alJ (119951) . 


We app ly a mass-loss enhancem ent for stars near critical ro¬ 
tation as in I Yoon & Langeil (120051) which considers a reduction 
of the critical rotational velocity for stars near their Eddington 
limit. It is still uncle ar whether rapid rotati on per se leads to an 
enhanced mass loss dMiiller & VinkI 1^141) . but it appears rea¬ 
sonable to consider that the mass-loss rate increases close to the 
Eddi ngton limit, wh ich is indeed reached sooner for rotating ob¬ 
jects (lLangeilll997l) . 


2.4. Model grid 



^ZAMS [^0] 


Fig. 1. Initial equatorial rotational velocity versus initial mass. Each dot 
in this diagram represents the evolutionary sequence in our model grid 
with the corresp onding initial param eters. Grey dots correspond to mod¬ 
els presented in iBrott et all (1201 lal) as well as previously unpublished 
models calculated by I. Brott, while black dots represent the 110 newly 
computed evolutionary sequences. 


3. Results 

3.1. Evolutionary tracks in the HR-diagram 

A selection of stellar evolution tracks is presented in Fig.|2] The 
luminosity of the stellar models for a given initial mass and ini¬ 
tial surface rotational velocity is shown as a function of their 
effective temperature. Tracks are shown for nine different initial 
masses from 60 M© to 500 M©, with initial surface rotational ve¬ 
locities of Okm/s, 400km/s, and 500km/s. 



Fig. 3. Evolutionary tracks of stars with an initial mass of 80 Mq, for 
initial rotational velocities of 0, 50, 100, 150, 200, 250, 300, 350, 400, 
450, and 500 km/s. Tracks with initial velocities that are even multiples 
of 50 km/s are drawn in blue, odd multiples in red. 


Figure [T] gives an overview of our grid of evolution models by 
indicating the initial masses and initial surface rotational veloci¬ 
ties of all computed model sequences. Because of the increase in 
the Eddington factor with mass, we decrease the maximum ini¬ 
tial rotational velocity for higher masses in order to avoid strong 
rotation ally induced m ass loss already on the zero-age main se¬ 
quence (lLangedl998l) . Whereas most sequences are computed to 
core-hydrogen exhaustion, some of the most massive and most 
rapidly rotating models were stopped shortly before their prox¬ 
imity to the Eddington limit caused numerical difficulties. 


Up to a luminosity of log L/L© - 6.5, at about 190 M©, the 
effective temperature of the ZAMS increases with increasing lu¬ 
minosity. For higher initial masses this behaviour changes. The 
ZAMS moves towards lower effective temperature with increas¬ 
ing luminosity as a result of stars having significantly increased 
radii and as a consequence of their proximity to the Eddington 
limit. This effect is discussed in more detail in Sect. 13.61 

Figure |2] shows only the tracks of the slowest and fastest ro¬ 
tators in our grid. However, we point out that in the HR-diagram, 
the evolutionary tracks do not change much below an initial ro- 
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Teff [kK] 


Fig. 2. Evolutionary tracks of massive stars during their core hydrogen burning evolution in the Hertzsprung-Russell diagram. For each selected 
initial mass (as labelled), tracks are shown for three different initial surface rotational velocities, uzams = 0,400,500km/s, in black, blue, and red, 
respectively. The time difference between two successive dots on each track is 10^ yr. The zero-age main sequence is drawn as a green dashed line. 
The end of the tracks corresponds to the terminal age main sequence. The approximate location of the Humphreys-Davidson limit is indicated by 
the black dashed line jHumphrevs & Davidsonlfr994l) . 


tational velocity of ~ 250km/s. This is demonstrated by the 
example of 8 OM 0 models in Fig. |3] for which only the tracks 
with initial rotational velocities of ~ 300km/s and higher de¬ 
viate significantly from the track of the model without rotation. 
As shown in Sect. 13.31 this is due to the absence of significant 
rotationally induced mixing of helium below the threshold ro- 
tati onal velocity for q uasi-chemically homogeneous evolution 
(cf. lBrott et al.l d^l lah ). We note that although the evolutionary 
tracks of the sequences computed with initial rotational veloc¬ 
ities below ~ 250km/s are almost identical, the corresponding 
models show significant differences concerning the evolution of 
the surface abundances of trace elements, e.g. boron and nitro¬ 
gen (cf. Sect. 13.3b . 

The models without significant mixing of helium, in the mass 
range between 6QMq and 8 OM 0 , expand during core hydrogen 
burning to surface temperatures as low as 5000 K. This occurs 
partly because of the relatively large amount of convective core 
overshooting in our models. A second reason is that when these 
models evolve through core hydrogen burning, they approach 
the Eddington limit because their L/M ratio is increasing, and at 
the same time their envelope opacity is increasing as cooler sur¬ 
face temperatures are achieved (cf. Sect. I3.61 l. For higher initial 
masses, the redward evolution is truncated at Teff ^ 2 5 000 Kbe- 
cause of the bi-stability mass-loss enhancement in the I Vink et al.l 


(l 200 (lh mass-loss recipe and the assumed mass-loss enhancement 
for stars near their Eddington limit (Sect.|3. 

Eigure|2] shows that the evolutionary tracks of even our 
slowly rotating models avoid the upper-right corner of the HR 
diagram. This is remarkable, since our very massive star models 
evolve very close to their Eddington limit, which leads to an in¬ 
flation of the envelope (cf. Sect. l3.61 l. We find that the high mass- 
loss rates at temperatures below ~ 30 000 K lead to significant 
helium enrichments for all stars above ~ 6QMq (cf. Sect. 13.3b 
such that, as core hydrogen burning continues, they evolve to¬ 
wards hotter rather than cooler surface temperatures. 

Eigure|2] also contains the empirical upper luminos- 
ity boundary of stars in the Milky Way, as derived by 
[Humphreys & DavidsonI (Il994h . As we see, our slowly rotat¬ 
ing models do penetrate the Humphreys-Davidson (HD) limit 
and spend a significant amount of time at cooler temperatures. 
Whether this prediction is in contradiction with observations 
for the EMC is currently unclear. The stellar statistics near the 
HD limit for EMC stars from published work is not very good 
(lEitzpatrick & Garmanvll9^ . And even in the Milky Way, stars 
above the HD limit are observed, the most prominent exam¬ 
ple being ? 7 Cari nae with a l uminosity of log L/Lq - 6.7 and 
Teff ^ 30 000 K dSmitSIMll)- 


Article number, page 4 ofl21l 







































K. Kohler et al.: Evolution of rotating very massive stars 


In any case, our models predict a short lived (~ 10^ yr) yel¬ 
low or red supergiant phase of stars in the mass range 60Mo 
to 80Mo during which the core is still burning hydrogen. Our 
models obtain a stellar wind mass-loss rate of the order of 
10 '^Moyr^^ during this stage. However, in this part of the HR 
diagram, mass-loss rates are very uncertain. We note that in par¬ 
ticular higher mass-loss rates (which may be due to the pulsa- 
tional instability of these models; Sanyal et al., in prep.) would 
lead to shorter life times in this evolutionary stage. 

The higher the initial mass of our model, the higher is their 
convective core mass fraction. While it exceeds 80% in zero-age 
main sequence stars above lOOM©, it reaches 90% at 3QQMq, 
and SOOMq stars with a convective core mass fraction of 95% are 
almost fully convective. This puts the convective core boundary 
far out into the stellar envelope, where the pressure scale height 
is rather small. Consequently, convective core overshooting does 
not have the same importance as at smaller stellar mass and plays 
only a minor role in most of the models presented here. 

The effect of rotation on the evolution of massive stars is 


discussed previously, for example in jMaeder&^leynei (2000^ 
[ le ger & Langed (l200(lh : lBrott et alJ (1201 lah : Chieffi & Limongil 
( 2013h . In our models, there are two main changes in the stel¬ 
lar evolution tracks in the HR diagram. First, the effective grav¬ 
ity is reduced as a result of the centrifugal acceleration, which 
leads to a decrease in the effective temperature and luminos¬ 
ity of a star compared to a non-rotating model. Second, above 
a threshold rotational velocity for which the timescale of rota¬ 
tional mixing becomes smaller than the nucle ar timescale, the 
stars evolve quasi-chemically homogeneously (lYoon & Langej 
l2005t IWooslev & Hegeiil2006ll . and the corresponding models 
evolve di rectly towards the helium main sequence in the HR- 
diagram dBrottet alJl2011ah . 

Both effects are clearly visible in our models. For stellar evo¬ 
lution tracks with initial rotational velocities below 250km/s, the 
first effect mentioned is dominant. It can be recognized most eas¬ 
ily from the ZAMS position of our stellar models (Figs.|2]and[3ll, 
where the models are dimmer and cooler the faster they rotate. 

The most rapidly rotating stellar models undergo quasi- 
chemically homogeneous evolution. As displayed in Fig. for 
t^zAMS = 400,500km/s, those with initial masses below ~ 
125Mo show a strong increase in luminosity and effective tem¬ 
perature. For more massive models, however, homogeneous evo¬ 
lution also leads to higher luminosities, but the surface tem¬ 
perature is decreased. Towards core helium exhaustion, strong 
Wolf-Rayet mass loss leads these stars to lower luminosities and 
larger surface temperatures (cf. Sect. I3.21 i. We return to the dis¬ 
cussion of quasi-chemically homogeneous evolution in Sect. l3.4l 
Isochrones in the HR-diagram based on tracks discussed here are 
presented in Appendix A. 


3.2. Mass loss and surface rotational velocity 

The evolution of the models for the most massive stars is 
strongly affected by mass loss. According to our mass-loss pre¬ 
scription (Sect. 13, our models undergo mass loss of t hree differ¬ 
ent strengths . Initially, the mass-loss rate proposed bv IVink et akl 
(l2000i 1200 ih is used. For our adopted composition, the mass- 
loss rate is initially of the order of 3 ■ IO^^Mq yr“' for the lOOM© 
models, while for the 500Mq models it is ~ 6 ■ lO^^Meyr^*. At 
this rate, only 10-20% of the initial mass will be lost over the 
lifetime of the stars (Fig.|4] right panel). 

The models that do not experience quasi-homogeneous evo¬ 
lution undergo an increase in their mass-loss rate at effective 
temperatures of ~ 25 000 K, which is referred to as the bi¬ 


stability jump (IVink et al.lll999l) . according to the Vink et al. 
prescription (Fig.|4]i- The corresponding mass loss is so intense 
that the models above IOOMq stop evolving redward at this 
stage, since their surfaces become strongly helium-enriched. For 
our most massive models, the 5OOM0 sequences, the maximum 
mass-loss rate at this stage is about 2 ■ IO^^Mq yr“' . 

Once the surface helium mass fraction has increased owing 
to mass loss to more than 40%, the Wolf-Rayet mass-loss rate 
is phased in, reaching its full strength at a surface helium mass 
fraction of 70%. As a consequence, our 100 - 5OOM0 models 
lose 50-80% of their initial mass before core hydrogen exhaus¬ 
tion. During the phase of Wolf-Rayet winds, the models evolve 
at decreasing luminosity and all accumulate in a narrow region 
in the HR-diagram (Fig.|3- 

The three different phases of mass loss can be clearly seen 
in both panels of Fig. |4] Again, for our most extreme models, 
the slowly rotating stars with an initial mass of ~ 500Mo, the 
largest obtained mass-loss rate is ~ 5 • lO“^M0yr“*, according 
to our Wolf-Rayet mass-loss prescription. At this moment, the 
stars have a luminosity of log L/Lq ^7.1. Assuming a terminal 
wind speed of 1000 km/s, we compute a wind-momentum-to- 
photon-momentum ratio of 77 ^ 2 for this situation. With the 
same numbers, we obtain a wind-kinetic-energy-to-luminosity 
ratio of 0. 003. The correspondin g wind darkening is therefore 
negligible dHeger & Langerl 19961) and significantly smaller than 
in Galactic Wolf-Rayet stars ( Lucv & Abbotd[r993h . 

The models that evolve quasi-homogeneously first increase 
their luminosity more than the models described above. Con¬ 
sequently, their mass-loss rate becomes larger than that of the 
inhomogeneous models. The homogeneous models do not reach 
the bi-stability limit as the increased helium surface abundance 
due to rotational mixing keeps their surface temperature above 
25 000 K. However, they can reach the Wolf-Rayet stage earlier, 
and in the end lose similar amounts of mass as the inhomoge¬ 
neous models. During the Wolf-Rayet stage, there are no clear 
characteristics that can distinguish between the two types of evo¬ 
lution. 



Fig. 5. Mass as a function of effective temperature for the models with 
initial masses of 60Mq, 70Mq, and SOMq, and for initial rotational ve¬ 
locities of 100, 400, and 500 km/s, as indicated. The time difference 
between two dots along each of the tracks is 10^ yr. 

Figure |5]compares the mass evolution of the homogeneously 
and inhomogeneously evolving 6QMq to 8OM0 models. It shows 
that in this mass range, the homogeneously evolving models lose 
more mass, with more than half of the mass being lost during the 
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Fig. 4. Mass as a function of effective temperature (left panel) and as a function of time (right panel). Stellar models that evolve towards effective 
temperatures below 27 000 K are affected by the bi-stability jump leading to an increase in the mass-loss rate. This is followed by a change in the 
slope of the mass as a function of time. The increase in mass-loss rate for rapidly rotating models is related to the change to the lHamann et al.l 
lll99.‘5h mass-loss recipe. 


Wolf-Rayet phase. The inhomogeneously evolving models lose 
almost equal amounts during their evolution through the hot part 
of the HR-diagram as at cool surface temperatures, where the 
latter mass loss occurs on a timescale of only 10^ yr. 

As a consequence of the initially rather moderate stellar wind 
mass loss, most of our models evolve at a constant rotational ve¬ 
locity for most of the core hydrogen burning. In th is case the loss 
of angular momentum through the stellar wind (lLangeilll998h 
and the increased momentum of inertia due to expansion is com¬ 
pensated by the transport of angular momentum from the con¬ 
tracting core to the envelope. As shown in Fig.|6] this holds even 
for the most massive stars in our grid. Only the models which un¬ 
dergo chemically homogeneous evolution, as a result of their in¬ 
creased luminosity and mass-loss rate, show a moderate decline 
of their surface rotational velocity before entering the phase of 
Wolf-Rayet winds. 



Fig. 6. The surface rotational velocity is shown as a function of time for 
stellar model sequences with four different initial masses (as indicated), 
and for three initial surface rotational velocities. The rapid decrease in 
rotation rate seen for c; 50,200 km/s is caused hy the enhanced mass 
loss at the hi-stability jump, whereas the steep decline of the rotation 
velocity for the fast rotators (cj =: 500 km/s) is a consequence of Wolf- 
Rayet type mass loss. 

The strong mass loss for Tsg < 25 000 K for the inhomo¬ 
geneously evolving models, and the Wolf-Rayet type mass loss 


for the homogeneously evolving ones both lead to a strong spin- 
down of the stars towards the end of their core-hydrogen burn¬ 
ing evolution. In fact, this i s only avoided for stars initially be¬ 
low ~ 30Mq (cf. Fig. 3 of I Vink et al.l (l2010l) ). Since the mag¬ 
netic coupling in our models ensures close-to-rigid rotation dur¬ 
ing core hydrogen burning, this implies that all our LMC mod¬ 
els above ~ SOMq lose so much angular momentum that they 
cannot be considered to produce candidates for long-duration 
gamma-ray bursts. In that respect, models that undergo chemi¬ 
cally homogeneous evolution of lower mass are better suited (cf 
Sect.[34li. 

Effective temperatures and radii of observed massive stars 
may be affected by optically thick stellar wind s, making the s tar 
appear larger and cooler. We use Eq. (14) from lLangeij (ll98^ 


r(R) = 


k\M\ 

AnR (v 

OO - Vo) Vo 


( 1 ) 


to estimate the optical depth of the stellar winds of our stel¬ 
lar models. Here, R designates the radius of the stellar model 
without taking the wind into account. This equation is derived 
from a y6-velocity law with (3 — 1. In this case, we use the elec¬ 
tron scattering opacity {k - cr(l -H X), where cr is the Thomson 
scattering cross section), an expansion velocity vq - 20km/s at 
the surface of the stellar model, and a terminal wind velocity of 
Voo = 2000 km/s. 

Eigure|7]shows the estimated optical depth of stellar winds as 
a function of time for several model sequences. The behaviour 
of the optical depths seen in the figure is mostly related to the 
change in the mass-loss rate. The optical depth increases for 
higher initial mass and higher rotation rate as a result of a corre¬ 
sponding increase in the mass-loss rate. While these numbers are 
only approximations, it shows that the winds of the most massive 
stars might already be optically thick (r > 1) on the zero-age 
main sequence, which implies that spectroscopically, these stars 
may already show Wolf-Rayet c h aracteristics at t his tim e (see 
also iGrafener & Hamamil (l2008h : ICrowther et al.l (120101) ). Eur- 
thermore, while the stars of WOMq and below are expected to 
show optically thin winds for most of the core hydrogen burning 
evolution, the optical depths of the winds for models for which 
a Wolf-Rayet type wind has been assumed may be quite large. 
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Fig. 7. Stellar wind optical depth according to Eq. (1), for some of our 
model sequences, shown as a function of time. We show eight stellar 
models with four different initial masses (60, 100, 300, 500 Afg), and 
two initial surface rotational velocities (0, 500 km/s) in black and red, 
respectively. The line for unit wind optical depth is plotted to facilitate 
the comparison. 


3.3. Surface abundance evolution 

Deriving masses or ages of stars using stellar evolution tracks or 
isochrones in the HR diagram is no longer straightforward when 
rotational mixing is important. A given pair of effective temper¬ 
ature and luminosity can be explained by several combinations 
of masses, rotational velocities, and ages. To uniquely determine 
the initial parameters and the age of a star, it is necessary to in¬ 
clude additional observables. The mass fractions or abundances 
of elements at the stellar surface are tracers of the rotational mix¬ 
ing. They can be used to uniquely identify the initial parameters 
and the age of a star in a three (or m ore) dimensional space of 
observables (cf. iKohler et al.l d2012h f. We therefore discuss the 
influence of time, initial rotational velocity and initial mass on 
the surface composition of our models. 

The rotational mixing in our models is sensitive to the gra¬ 
dient of the mean molecular weight ji. Essentially, any signifi¬ 
cant /i-barrier prevents rotational mixing. Therefore, mixing of 
helium can only occur in our models as long as they are quasi- 
chemically homogeneous. This, in turn, requires the timescale 
for rotational mixing to be shorter than the nuclear timescale, as 
the star attempts to establish a p-barrier on the nuclear timescale. 

In the fastest rotators, e.g. in the IOOMq and BOOM© models 
initially rotating with ~ 400 and 500 km/s, the surface helium 
abundance goes almost all the way to Tj = 1 towards core hydro¬ 
gen exhaustion (cf. Fig.[^. However, the same figure shows that 
this evolution is truncated for the rapidly rotating 6QMq mod¬ 
els, at Ts = 0.68 and - 0.80 for initial rotational velocities 
of ~ 400 and 500 km/s, respectively. These models spin down 
during core hydrogen burning (Fig. |6]l such that the timescale 
for rotational mixing eventually becomes larger than the nuclear 
timescale. 

Figure [8] shows that the non-rotating SOOMq sequence also 
evolves to Ts = 1. The reason is the large convective core fraction 
of this model, and its large mass-loss rate. Whereas the model 
always keeps a small radiative envelope mass, the equivalent 
amount of mass is lost on a short timescale such that core and 
surface abundances become almost equal (cf. Sect. 13.41 Eq. 2). 
For the rapidly rotating IOOMq and SOOMq models, mass loss 
rather than rotational mixing must take over to enforce chemical 


homogeneity, since these models also spin down quite dramati¬ 
cally (Fig.|3. 

Figure|3(left panel) demonstrates the role of rotational mix¬ 
ing of helium in greater detail using the example of 80Mo se¬ 
quences. In the models initially rotating with 0, 50,100,150, and 
200 km/s, such mixing is essentially negligible, and only mass 
loss can increase the helium surface abundance close to core 
hydrogen exhaustion, which is why all of these models evolve 
along the same track in the HR-diagram (Fig. (3). The models be¬ 
tween 250 and 400 km/s start out with chemically homogeneous 
evolution (though the one with Umt.i = 250 km/s only for a short 
amount of time), and truncate the homogeneous evolution later 
in time for higher initial rotation. The models which initially ro¬ 
tate with about 450 and 500 km/s undergo chemically homoge¬ 
neous evolution all the way (though with the help of mass loss in 
the end), and their evolution in the HR diagram is also practically 
identical (Fig. [3. 

The evolution of the surface nitrogen abundance in our mod¬ 
els of rotating and mass losing stars follows different rules. Dur¬ 
ing a phase of chemically homogeneous evolution, nitrogen is 
quickly mixed from the core to the surface establishing the CNO- 
equilibrium value in atmospheric layers (Figs. 0 and 0 left pan¬ 
els). Mass loss can also lead to nitrogen enhancements, even 
without rotational mixing (Fig.[3- Figure |3 shows that a strong 
nitrogen surface enhancement can be obtained even in models 
with initial rotational velocities well below the threshold value 
required for chemically homogeneous evolution. For example, 
the 80Mo models with iifot.i = 150 and 200 km/s enrich nitrogen 
at the stellar surface by factors of 4 and 12 by rotation, i.e. be¬ 
fore mass loss kicks in. This order of magnitude of nitrogen en¬ 
richment is as expected, since our models were calibrated to in¬ 
creases the surface nitrogen abundance by a factor of about 3 for 
stars of 13-15M0 for an initial rotational velocity of 150km/s, 
and rotational mixing is stronger in more massive stars. 

The reason is that nitrogen, carbon, and oxygen, with a max¬ 
imum mass fraction of less than one per cent, are just trace ele¬ 
ments, such that even significant internal gradients do not lead to 
strong p-barriers. Therefore, at the beginning of core hydrogen 
burning, before significant amounts of hydrogen have been con¬ 
verted to helium, but after the CNO-cycle has already strongly 
enhanced the nitrogen abundance in the convective core, rota¬ 
tional mixing can bring nitrogen out of the core into the radiative 
envelope and later on also to the surface. 

The surface enrichment of helium and nitrogen explained 
above allows a simple understanding of the occurrence of sur¬ 
face enrichments in the HR diagram. As shown in Fig.[T0l stars 
initially rotating with rotational velocities of 200 km/s or less 
that are less luminous than log L/Lq ^ 6.2 are not expected to 
show any helium surface enrichment (i.e. < 0.28) as long as 

their surface temperature is higher than ~ 20 000 K. The same 
is true for stars initially rotating slower than 300 km/s below 
log LJLq ^ 5.6, and for those initially rotating slower than 
400km/s below log L/Tq - 4.8. Similar higher thresholds can 
be read off from Fig.fTOlfor larger surface helium mass fractions. 
Given that iRamfrez-Agudelo et al.l (1201 3h found that 75% of all 
O stars rotate slower than 200 km/s and that amongst the 31 02 
to 05 stars none was found to rotate faster than ~ 300 km/s, the 
helium enrichment in EMC stars below ~ IOOMq is expected to 
be quite small during the main-sequence phase. 

In contrast, as more than half of all O stars were found to ro¬ 
tate faster than 100 km/s, nitrogen enrichment by at least a factor 
of 2 is expected to be almost ubiquitous above log L/Tq - 6.0, 
and quite frequent in main-sequence O stars in general (Fig.fTTTi. 
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Fig. 8. Helium (left panel) and nitrogen (right panel) surface mass fraction as a function of time for models of 60, 100, and 300 Mq and rotation 
rates explained in the figure key. While for the fast rotators the enhancements are mostly due to rotational mixing, for the slowly rotating models 
the increase in helium and nitrogen in the 100 Mq and 300 Mq mass models is solely due to mass loss. 




Fig. 9. Helium (left panel) and nitrogen (right panel) surface mass fraction as a function of time for models of 80 Mq with approximate initial 
rotational velocities of 0, 50, 100, 150, 200, 250, 300, 350, 400, 450, and 500 km/s, during their core hydrogen burning evolution. 


Changes in the amount of neon, sodium, aluminum, and 
magnesium at the surface also occur, indicating that in very mas¬ 
sive stars the MgAl- and the NeNa-cycles are active. In our mod¬ 
els, the amount of magnesium is reduced while aluminum is pro¬ 
duced. The mass fraction of neon increases in our models, while 
sodium decreases. For more details, we refer to the electronic 
data published with this paper. 


Lithium, beryllium, and boron are destroyed in the stel¬ 
lar interior (IMc William & Rauchl 120041) at temperatures above 
~ 2.5 • 10® K for lithium, ~ 3.5 ■ 10® K for beryllium, and 
~ 5.0 ■ 10® K for boron. In our models, temperatures below 
5 ■ 10® K are only found in the outer envelope. Therefore lithium, 
beryllium, and boron can exist in the outer envelope, while they 
are destroyed in deeper layers. For the slowest rotating models, 
the surface abundances of these three elements remain constant 
until layers which were exposed to higher temperatures earlier in 
the evolution are exposed to the surface due to mass loss. When 
that happens, the mass fractions of beryllium, boron, and lithium 
are quickly reduced. For fast rotating models, rotational mixing 
and mass loss lead to a gradual decrease in lithium, beryllium, 
and boron over time. 



t [Myr] 


Fig. 12. Kippenhahn-diagram for our 200 Mq model with Vzams - 
300 km/s. The black solid line gives the stellar mass as a function 
of time. Blue shading indicates thermonuclear energy generation (see 
colour bar to the right of the plot). Green hatched parts show convective 
regions, and convective core overshooting is indicated in red. Three dif¬ 
ferent regimes can be distinguished according to the rate at which the 
mass of the convective core decreases (see text). 
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Fig. 10. Evolutionary tracks in the HR diagram of stellar models initially rotating with approximately 200 km/s, and with initial masses of 12, 15, 
20, 25, 30, 40, 50, 60, 70, 80, 100, 125, 150, 175, 200, 230, 260, 300, 400, and 500Mq. Overlaid are lines (in blue) of constant helium surface 
mass fraction for Y=0.28, 0.35, 0.50, and 0.70. Lines of constant helium surface mass fraction for Y=0.28, 0.35 corresponding to models with 
approximate initial rotational velocities of 300 and 400 km/s are also shown. 



Fig. 11. Evolutionary tracks in the HR diagram of stellar models initially rotating with approximately 200 km/s, and with initial masses of 6, 7, 
8, 9, 10, 12, 15, 20, 25, 30, 40, 50, 60, 70, 80, 100, 125, 150, 175, 200, 230, 260, 300, 400, and 500Mo. Overlaid are lines (in blue) of constant 
nitrogen surface mass fraction, corresponding to nitrogen enhancement factors of 2, 5, 10, and 20, as indicated in blue. Lines of constant nitrogen 
enhancement factors of 2 and 3 corresponding to models with approximate initial rotational velocities of 150 and 100 km/s are also shown. 


3.4. Chemically homogeneous evolution 

As mentioned in Sect. 13.11 we can divide our models into three 
classes. The first one (which we call Class O) likely corresponds 
to most stars in observed stellar samples and contains the models 
that evolve in the normal way, i.e. in which the rotationally in- 
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Fig. 13. Central Yc and surface Kj helium abundance of our 200 Mq 
model with Vzams - 300km/s (cf. Fig. [13, and the ratio Ks/Fc as a 
function of time. 


duced mixing of helium is negligible. The second one describes 
the models that undergo quasi-chemically homogeneous evolu¬ 
tion (Class H), which correspond to the initially fastest rotators. 
As the third class, we have the models that start out evolving ho¬ 
mogeneously, but which spin down such that the rotational mix¬ 
ing of helium stops (Class HO), after which they evolve redward 
in the HR diagram as the ordinary models. These models re¬ 
tain a memory of their past homogeneous evolution, in that they 
will keep an enhanced he lium surface a bundance and a higher 
luminosity-to-mass ratio (lLangeiiri992h compared to ordinary 
stars throughout the rest of their core hydrogen burning evolu¬ 
tion. 

At the highest masses considered here {M >15QMq), our 
models may also undergo quasi-chemically homogeneous evo¬ 
lution without rotationally induced mixing, which is due to a 
combination of an extremely high fraction of the convective core 
to the total stellar mass and a very high Wolf-Rayet type mass- 
loss rate. We illustrate this by the example of our 200 Mq se¬ 
quence with I'zAMS - 300km/s in Fig. For such very massive 
stars, the convective core comprises a major fraction of the stel¬ 
lar mass. 

We can divide the core hydrogen burning evolution of this 
model into three parts according to the mass evolution of the 
convective core. During the first part, the mass of the convec¬ 
tive core decreases as a function of time, but more slowly than 
the mass of the star, such that the mass of the radiative enve¬ 
lope decreases and the convective core mass fraction increases. 
During this phase, the model undergoes quasi-chemically homo¬ 
geneous evolution (cf. Fig.lOll. At f ^ 1.2 Myr, the model tran¬ 
sitions to ordinary evolution, from which time on the convective 
core mass decreases somewhat faster than the total mass. During 
this second part of its evolution, the helium surface abundance 
increases, even though at a much lower rate than the central he¬ 
lium abundance (Fig. [Ell. Finally, in part three, the Wolf-Rayet 
mass loss kicks in, which leads to a very small mass of the non- 
convective stellar envelope (Menv - 5Mq). The amount of mass 
corresponding to the non-convective envelope is lost (at a rate 
of ~ 2 10^'^Moyr^') on a timescale of ~20000yr which corre¬ 
sponds to only about 1% of the core hydrogen burning time, th. 
Consequently, the surface helium mass fraction during this stage 
is roughly equal to (within 1%) the central helium mass fraction. 


Thus, regardless of rotation, the model is extremely chemically 
homogeneous during this phase. 

We define a fourth type of evolution (Type M) by 


iffgnv 

M 


< 0.1 Th, 


( 2 ) 


where th is the hydrogen burning timescale. Since this condi¬ 
tion ensures that the surface-to-central helium abundance ratio 
is above ~ 0.9 — i.e. similar to or smaller than in the case of 
rotationally induced quasi-chemically homogeneous evolution 
— we can subdivide our models into further classes. For ex¬ 
ample, the 200 Mq sequence discussed here can be considered 
of Class HOM, as it first undergoes chemically homogeneous 
evolution due to rotational mixing, then ordinarily, and finally 
chemically homogeneous evolution due to mass loss. 

It is instructive to consider Fig. [T4lto understand which evo¬ 
lutionary classes are realized by our models. The non-rotating 
sequences with initial masses of 300Mq and 500Mq achieve 
chemical homogeneity at a central helium mass fraction of 0.85 
and 0.75, respectively, and thus belong to Class OM. The non¬ 
rotating IOOM0 sequence does not reach homogeneity. Conse¬ 
quently, although it achieves a quite high final surface helium 
mass fraction of ~0.65, it belongs to Class O. The most rapidly 
rotating models depicted in Fig.[T4lnever show a significant dis¬ 
crepancy between their central and surface helium abundances, 
which implies that they evolve from H-type to M-type evolu¬ 
tion and belong to Class HM. The stars rotating initially with 
~ 300km/s all undergo an H—>0 transition, where the two more 
massive ones move further on to Type M evolution. As stars in 
the Classes HOM and HM both start out and end their chemi¬ 
cally homogeneous evolution, we will not distinguish them fur¬ 
ther and designate them both as Class HM. 



Fig. 14. The mass fraction of helium in the stellar core Y^ and at the 
surface Y^ are depicted for stellar models of 100, 300, and 500 Mq and 
rotation rates explained in the figure key. For a given rotation rate, ini¬ 
tially more massive stars increase their surface helium abundance more 
quickly. Every 10^ yr the models are highlighted by filled circles. Chem¬ 
ically homogeneous evolution is indicated by equal changes in both 
mass fractions and corresponds to a slope unity. At Y^ > 0.8, all cal¬ 
culated stellar evolution models show an evolution toward the line of 
equal Y^ and Fj, caused by mass loss. 

Thus, we have five different classes of evolution, H, HM, 
HO, OM, and O, where all but the last involve quasi-chemically 
homogeneous evolution. Class O can be identified in Fig. [TSl 
(right panel), which shows the close-to-final surface helium mass 
fraction for all our model sequences. Looking at the slowly ro¬ 
tating models, this figure shows that mass loss starts affecting 
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Fig. 15. The grid of all i nitial masses M/ams and surface rotational velocities Vzams is shown, including the stellar models with initial masses 
above 19 Mq published in iBrott et al.l ll201 lah . The colour coding corresponds to the surface rotational velocity (left panel) and the helium mass 
fraction at the surface (right panel) at the time when the central helium mass fraction has reached 95%. 


the surface helium abundance above an initial mass of ~ 65Mq 
(log MIMq ^ 1.8), and that M-type chemically homogeneous 
evolution is obtained above ~ I6OM0 (logM/M© ^ 2.2). Ac¬ 
cording to our definition, the slow rotators below ~ I6OM0 
are thus in Class O, the more massive ones in Class OM. For 
the fast rotators, the dividing line between homogeneous and 
inhomogeneous evolution depends on stellar mass. At 2QMq 
(logM/M0 ^ 1.3), stars rotating initially faster than ~ 500 km/s 
evolve homogeneously (Class H), while at ~ I6OM0 the crit¬ 
ical velocity is at ~ 350 km/s, and more massive stars above 
this velocity are in Class HM. Below the dividing line defin¬ 
ing Class H, between 2QMq and I6OM0 is a stretch of Class HO 
models, which comprises a larger initial mass range for higher 
initial mass. 

As mentioned in Sect. 13.11 rapidly rotating stellar models 
with initial masses greater than 125 M0 show significantly dif¬ 
ferent behaviour in the HR diagram (see Fig. |2) to the stellar 
models with M < 125 Mo (see also Sect. I3.6l l. Nevertheless, 
they behave in a similar way in the Kippenhahn diagram and 
when comparing the helium mass fraction in the core and at the 
surface. 

Stellar evolution models above 80 Mq have strong mass loss. 
Ev en the fastest rotators are slowed down significantly when 
the iHamann et al.l (Il995h mass-loss rate is applied. For models 
above Mzams - 100 Mo, the rotation rate at 95% helium mass 
fraction in the core is below 100 km/s. The surface rotational ve¬ 
locity goes down to less than 50 km/s for models above 150 Mo, 
independent of the initial surface rotational velocity (Fig. [Bl 
right panel). 

From the above, we can come to several conclusions which 
are relevant in comparison to observed stars. First, helium- 
enriched single stars below ~ 65Mq are prime suspects of H- or 
HO-type chemically homogeneous evolution, whereas M-type 
evolution can be excluded for them. These stars are expected 
to preserve their rapid rotation throughout their core hydrogen 
burning evolution. As a consequence, we may expect some cor¬ 
relation of the helium surface abundance with the stellar rotation 
rate in the considered mass regime. 

For very massive stars (above ~ 2QQMq) rotation, whether 
initially fast or not, makes little difference, i.e. the evolution 
becomes almost independent of the initial rotational velocity. 
These stars are efficiently spun down by mass loss, and all un¬ 


dergo chemically homogeneous evolution during their advanced 
hydrogen-burning evolution. Their surface helium abundance is 
always very close to their central helium abundance because of a 
rapid loss of the thin outer radiative envelope through mass loss. 
As a consequence, these sequences can produce models that have 
a very high surface helium mass fraction (Tj > 0.9) and are still 
undergoing core-hydrogen burning. 

3.5. Mass-luminosity relation 

Stars of higher mass are increasingly luminous. The most mas¬ 
sive models in our grids, at 500Mo, radiate at more than IO^Lq. 
Figure [16] shows the evolution of selected non-rotating and of 
rapidly rotating model sequences in the mass-luminosity plane, 
from core hydrogen ignition to core hydrogen exhaustion. At the 
lowest considered masses (15Mo), the models evolve vertically 
upward since their mass loss is negligible. The more massive, 
rapidly rotating models evolve to higher luminosity, and they 
turn left toward lower mass. 

A comparison with the mass-luminosity relatio ns for chemi¬ 
cally homogeneous stars of iGrafener et al.l (1201 Ih reveals that, 
except for the depicted 15Mo model, the rapidly rotating se¬ 
quences shown in Fig. [16] undergo quasi-chemically homoge¬ 
neous evolution. They start at the mass-luminosity (ML) relation 
for zero-age main sequence stars and finish close to the ML- 
relation for chemically homogeneous helium stars, indicating 
their nearly homogeneous chemical structure. Th is also demon¬ 
strates that the agreement of the ML-relations of iGrafener et al.l 
(1201 ih (which represent fits to stellar models computed with a 
different code) with our models is very good. 

We consider a mass-luminosity relation for chemically ho¬ 
mogeneous stars as L ~ where jj. is the mean molecular 

weight of the stellar gas. The power-law exponent a thus de¬ 
scribes the sl ope of the ML-relation in the lo g M - logL-plane. 
According to iKippenhahn & Weigertl (Il990h . we have a > 1, 
and O' ^ 1 for M — » 00 , due to the increasing radiation pres¬ 
sure for higher masses (Fgas/(Fgas + Frad) ^ 0). The inequality 
or > 1 implies that, in the frame of simple opacity laws, massive 
main sequence stars will never exceed the Eddington limit. The 
limit of a —> 1 implies that for sufficiently massive stars, their 
hydrogen burning life time th - EjE becomes independent of 
mass, since their nuclear energy reservoir E is proportional to 
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3.7, and yUtams//Uzams - 2.23), the corresponding factor is only 2.4 
at SOOMq. 


3.6. Near the Eddington limit 

In its general form, the Eddington limit is complex as it involves 
an appropriate mean of the total of all opacity sources. This total 
opacity k may be a function of depth in the stellar atmosphere; 
therefore, it is non-trivial to uniquely define the circumstances in 
which the star encounters its Eddington limit, nor to assess the 
consequences such an encounter may have. The proximity of a 
star of mass M and luminosity L to the Eddington limit is usually 
expressed in terms of the Eddington factor 


r = 


kL 

AncG M 


'^Zlk 

eg 


(3) 


Fig. 16. Luminosity as a function of mass for selected non-rotating 
(blue) and rapidly rotating models (Vmij = 500km/s; red). The labels 
indicate the initial mass of the considered sequences. The tracks end at 
a central helium mass fraction of 17 = 0.98. Overplotted are the mass- 
luminosity relations of iGrafener et all ll201 ih for chemically homoge¬ 
neous stars with a hydrogen mass fraction of X=0.74 (labelled ZAMS) 
and X=0 (labelled TAMS). The straight line labelled “L ~ M” indicates 
the smallest expected slope of the mass-luminosity relation. 



Fig. 17. Exponents of the mass-l uminosity relation of the form L ~ 
for homogeneous stars from IGrafener et ahl ll201 ih . for a chemi¬ 
cal composition corresponding to our zero-age main sequence models 
(labelled with the subscript ZAMS), and for pure helium stars of the 
corresponding metallicity (labelled with the subscript TAMS). Here, p 
is the mean molecular weight of the stellar gas, which is computed as 
l/p = 2X + 0.75K-t- 0.5Z, with X, Y, and Z being the hydrogen, helium, 
and metal mass fraction, and assuming complete ionization. 


their mass, just as their energy loss rate E - L. Figure [16] shows 
that our most massive models reach cr-values very close to one. 
We find the limiting stellar life time to be close to 1.9 Myr (cf. 
Table Bl). 

Figure [T7| also shows that, at a given mass, the slope of 
the mass-luminosity relation is smaller for the homogeneous 
hydrogen-free stars, compared to that for the zero-age main se¬ 
quence stars because the helium stars are much more luminous 
and thus more radiation pressure dominated. Furthermore, the 
power-law exponent fi also approaches one for the highest con¬ 
sidered masses. While a \5Mq helium star is about 20 times as 
luminous as a ISMq zero-age main sequence star (yS(15Mo) - 


where the constants have their usual meaning and is the ef¬ 
fective temperature and g the surface gravity. 

Considering only the dominant contributor to the opacity, i.e. 
photon scattering on free electrons, greatly simplifies the concept 
of Eddington limit as electron scattering is a grey process and, 
for the hot stars considered here, independent of depth in the 
atmosphere. For a star of mass M and luminosity L we obtain 
Eg = KsLKAttcGM), where A-g = crg(l + X) is the opacity due to 
Thomson scattering. It holds that T > Eg as a > /Cg. Even in the 
most massive stars, the Eddington factor Eg will not exceed unity 
(see Sect. [331 i. Even so, since the mass-luminosity exponent a > 
1, main sequence stars will get ever closer to the Eddington limit 
for electron scattering the higher their mass. 

Considering the true Eddington factor, we must expect that 
the situation E = 1 is actually achieved in very massive stars. 
Indeed, while Ke - 0.34 in our zero-age models, their true surface 
opacity is typically A'surface - 0.5. We can thus expect that E = 1 
is achieved at Eg ^ 0.7. 

In Fig. [T8| we plot the quantity ;= T'^f^lg - c/(A'gCr)Eg 
as a function of the effective temperature of selected evolution¬ 
ary sequences. We note that .if is normalized to the solar value, 
logjZ^ ^ 10.61 for convenience. For instance, for the initial 
composition of our stars we have log .if/.if© - 4.6 -H log Eg. 
The quantity .if is proportional to Eg and can be derived from 
spectroscopi c observation of stars with out knowledge of their 
distance (cf. iLanger & Kudritzkl (l2014lH . 

Since for stars of constant mass .if ~ L, the quantity .if 
behaves in a similar way to the stellar luminosity, and the evo¬ 
lutionary tracks in the .if-Tgff-diagram partly resemble those in 
the Hertzsprung-Russell diagram. However, unlike the latter, the 
.if-Tgff-diagram has an impenetrable upper limit; the Eddington 
limit. 

Comparing Fig.[T8]with Fig.|2| we see that while the tracks of 
the slow rotators in Fig.|2|are almost horizontal, the correspond¬ 
ing tracks in Fig. [18] are significantly steeper. Their luminosi¬ 
ties are nearly constant despite strong mass loss, because their 
increasing mean molecular weight p compensates for their de¬ 
creasing mass in the mass-luminosity relation (Sect. [33]| . While 
it remains hidden in Fig. ]2] Fig. [TSl reveals nicely that the mass 
loss drives these stars towards the Eddington limit. The same is 
true for the rapidly rotating stars; their tracks move steeply up¬ 
ward in both diagrams. 

In Fig. [18] we have drawn the horizontal line indicating 
Eg = 0.7 for our initial chemical composition — which corre¬ 
sponds to log.jSf/.5f^ ^ 4.445. From the discussion above, we 
should expect that stellar models with unchanged surface com¬ 
position should not be found above this line (as long as their 
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Fig. 18. The evolutionary tracks in the ^ — Teff-diagram, where ^ = T^glg, for models initially rotating with lOOkm/s (for 60, 80, 100, 125, 
150, 200, 300, and 500Mo; alternating red and blue lines) and with 400km/s (100, 150, 200, and SOOMg; green lines). The right y-axis shows 
the Eddington factor for electron scattering opacity T^, which is proportional to for a hydrogen mass fraction of X = 0.74. For X = 0, 
logT,, is larger by 0.23. Triangles, pentagons, and heptagons mark the locations where the surface helium mass fraction reaches 0.3, 0.5, and 0.7, 
respectively, for the presented evolutionary sequences. The dotted horizontal lines marks the value of Tg = 0.7 for X = 0.74 (green) and 2f = 0 
(blue), approximately identifying the true Eddington limit that cannot be exceeded (see text). 


hydrogen and helium can be considered fully ionized). This is 
confirmed by the triangles placed on the evolutionary tracks, in¬ 
dicating the position at which a surface helium mass fraction of 
Y - 0.3 is achieved, and which are all found below Tg = 0.7. 
Similarly, even the models with the most helium-enriched sur¬ 
faces stay below the upper horizontal line for re(2(' = 0) = 0.7; 
in other words, the diagram shows that it is not the electron¬ 
scattering Eddington limit which constrains the evolution of the 
most massive stars, but the true Eddington limit. Eurthermore, 
Eig.HH indicates that the limiting electron-scattering Eddington 
factor — which can be read off from Eig. [T8]to be a little bit be¬ 
low Eg = 0.7 — is independent of the hydrogen/helium surface 
abundances. 

We find our models near the Eddington limit to undergo a 
significant envel ope inflation . T his phenomenon h as been de¬ 
scribed before bv iKatol (Il986h and llshii et al.l (1 19991) f or zero-age 
main s e quence stars o f differ e nt metallicities, and byllshii et alJ 
(Il999l) : iPetrovic et al.l (l2006l) : iGrafener et alJ (120121) for Wolf- 
Rayet stars. Eigure [T^ shows the location of our zero-age main 
sequence stars, and of zero-age helium stars with the same metal- 
licity, in the HR-diagram. Both curves show a maximum effec¬ 
tive temperature, which is reached at about SOM© (140 000 K) 
for the helium stars, and at 200Mq (54 000 K) for the hydrogen- 
rich stars. Above these masses, the effective temperature de¬ 
creases again owing to the envelope inflation. 



logWK 

Fig. 19. The zero-age main sequence of our non-rotating stellar mod¬ 
els in the indicated mass range (solid red line), compared with zero-age 
main sequences helium star models for Milky Way (dashed blue line) 
and SMC (dash-dotted blue line) composition. Labelled dots along the 
lines imply masses in solar units. Owing to envelope inflation the helium 
main sequences bend towards cooler temperatures close to the Edding¬ 
ton limit. 


As seen in Eig. [191 the line describing the helium stars bends 
sharply to cooler temperatures at about 1 000 OOOL©, and is 
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even e xpected to cross the hydrogen-ZAMS (see also llshii et al.l 
(Il999l) ). implying that more luminous helium stars are cooler 
than equally bright zero-age main sequence stars. This has 
the consequence that stars undergoing chemically homogeneous 
evolution above ~ \25Mq evolve to cooler surface temperature 
(see Fig.|2]i. 

The envelope inflation occurs because layers inside the stel¬ 
lar envelope reach or exceed the true Eddington limit. The infla¬ 
tion begins before the critical value of Tg - 0.7 is reached, i.e. 
the electron-scattering Eddington factors are only Eg - 0.31 and 
Tg ^ 0.41 for the \5Mq helium star and \9QMq zero-age main 
sequence star, respectively, because of the opacity peaks of iron 
and helium which can reach values of 1 cm^/g or more. While 
the subsurface layers with such opacities always turn convec- 
tively unstable, convection, according to the mixing-length the¬ 
ory, is mostly very inefficient in transporting energy in these lay¬ 
ers. Consequently, the high radiation pressure pushes the over- 
lying layers outwards, until the opacity drops such that the true 
Eddington factor falls short of T = 1. 

Many of these models also develop densit y inversions nea r 
the layer with the maximum Eddington factor (IJoss et aOl973h . 
The corresponding inward directed gas pressure gradient then 
allows the star to retain layers whose Eddington factor exceed 
the critical value of E = 1. While the stability of these structures 
remains to be investigated, any instability which would tend to 
iron out the density inversion will likely lead to a further infla¬ 
tion of the overlying layers. A further degree of complexity is 
added by the fact that we And many models close to the Edding¬ 
ton limit to be unstable at least to radial pulsations (the only pul¬ 
sation mode which we can see with our one-dimensional hydro- 
dynamic stellar evolution code), and that convective velocities 
close to the sound speed m ay imply significant acoustic fluxes 
(iGoldreich & Kumar! 1990h . A deeper investigation of these phe¬ 
nomena exceeds the scope of the present paper and will be pur¬ 
sued elsewhere. 

4. Comparison with previous resuits 



Fig. 20. Compar ison of the ZAMS position of the models from 
I Yusof et ^ ll201.3li which include rotation (blue dots) with our tracks 
for Urot.i - 300km/s (red lines) in the HR diagram. We show tracks 
with initial masses of SOMq, IOOMq, 125Mo, 150Mo, 175Mo, 200Mq, 
230Mq, 260Mo, 3 OOM 0 , 400Mq, and 500Mq. The time interval between 
two dots on the tracks corresponds to 10^ yr. 

Recently. I Yusof et'aP (1201 3h published a grid of stellar evo¬ 
lution models for very massive stars. Eight of their evolution¬ 
ary sequences were modelled using a similar initial composition 



Fig. 21. Stellar mass at the end of core hydrogen burning versus ini¬ 
tial mass for our sequences with initial rotational velocities of Vrot.i - 
Okm/s, 300km/s, and 500km/s (filled symbols). T he open symbols cor¬ 
respond to the LMC models of lYusof et al.l (1201.30 with and without ro¬ 
tation. The dotted straight line corresponds to the location where the 
stellar mass at the end of core hydrogen burning equals the initial mass, 
i.e. to negligible mass loss. 

to the one we used here, except that Yusof et al. used a scaled 
solar metal mix, while we take the meas ured LMC composi¬ 
tion where possible (cf. lBrott et al.l(l2011ah '). Yusof et al. present 
non-rotating models with initial masses of I2OM0 , ISOMq, and 
5OOM0 as well as rotating models at I2OM0, ISOM©, 200Mq, 
3OOM0, and 5OOM0. The rotating models of Yusof et al. had ini¬ 
tial rotational velocities of 400km/s or more (cf. their Table 2). 

4.1. Subsurface convection and envelope inflation 

A major difference between the models of Yusof et al. and our 
models occurs because of the different treatment of convection in 
the stellar envelopes. Starting at zero age, all our models contain 
subsurface convection zones which occur because o f the opacity 
peaks of iron and helium (cf. ICantiello et al.l (l2009h l. When us¬ 
ing the standard mixing-length theory, where the mixing length 
is assumed to be proportional to the pressure scale height, we 
And the convective energy transport in these zones to be mostly 
inefficient because of the small local thermal timescales imposed 
by the low densities. As a consequence, we obtain a rich enve¬ 
lope phenomenology as described in Sect. I3.6I including enve¬ 
lope inflation and density inversions. 

Yusof et al. describe this situation as unphysical, and avoid 
it by assuming the convective mixing length to be proportional 
to the density scale height. As a consequence, when the star at¬ 
tempts to establish a density inversion, the density scale height 
tends to infinity, which imposes the convective energy transport 
to be extremely efficient, and opacity peaks will have few conse¬ 
quences. While the exi stence of inflated enve l opes in nature re¬ 
mains to be shown (cf. iGrafener et al.l (l2012h : iGrafener & VinkI 
dmH)), we consider it likely, since the short thermal timescales 
of the envelopes of very massive stars cannot be in accordance 
with efficient convective energy transport. Wfiile tfiis will surely 
be studied in more detail in the near future, we restrict ourselves 
here to pointing out the consequences of the different assump¬ 
tions on the convective efficiency. 

Strong differences between the models of Yusof et al. and 
our models already occur at the zero-age main sequence. As de¬ 
scribed in Sect. I3.1I our zero-age models (with our metallicity 
of Z = 0.0047) reach a maximum effective temperature of about 
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56 000 K at ~ IOOMq, which agrees very well with the results 
of llshii et alJ for Z = 0.004. In Fig. |20l we compare the 

location of the rotating zero-age models of Yusof et al. (accord¬ 
ing to their Table 2) with the early hydrogen burning evolution 
of our models with an initial rotational velocity of ~ 300km/s. 
Above ~ 120Mo, the models of Yusof et al. are hotter than ours, 
and the more so the larger the mass. 

This discrepancy is likely explained by the difference in the 
treatment of the envelope convection, since our treatment yields 
a larger radius inflation the closer the star is to its Eddington limit 
fSect.lra. At a given initial mass, the discrepancy between the 
models of Yusof et al. and our models is therefore expected to 
increase with time at least initially, since the models increase 
their luminosity and decrease their mass (cf. Fig. fTbl) . and thus 
their Eddington factors grow. The biggest consequence of this is 
that — although we use a very similar mass-loss prescription to 
that of Yusof et al. — our mass-loss rates are higher because our 
stars have cooler surfaces. 

4.2. Mzams versus Mjams relation 

Figure |2T] shows that as a result of the different treatment of en¬ 
velope convection in our models and those of Yusof et al., our 
non-rotating models above ~ ISOMq undergo a dramatic mass 
loss during core hydrogen burning, leading to masses at core hy¬ 
drogen exhaustion of the order of SOMq even up to the highest 
initial mass of SQQMq. Yusof et al. did not consider non-rotating 
models between ISOM© and SQQMq. However, their non-rotating 
models below ISOMq and at SOOMq remain close to the diago¬ 
nal in Fig.|2T] While our SQQMq model ends hydrogen burning 
with a mass of 87Mo, theirs does so with 2S9Mq. We discuss 
the consequences for the final fate of the stars in Sect. a Re¬ 
markably, because of the low mass-luminosity exponent of the 
most massive stars considered here, the core hydrogen burning 
life time of the SQQMq sequence of Yusof et al. (1.9 Myr) resem¬ 
bles ours, indicating that the life time of the most massive stars 
is not strongly affected by the differences discussed above. 

The non-rotating SQQMq sequence of Yusof et al. ends core 
hydrogen burning as a nearly chemically homogeneous helium 
star with a surface helium mass fraction of 0.97. Our Fig. [19] 
above shows that at the considered metallicity, our helium star 
models start to h ave i nflated envelopes abov e ~ SOM© (see also 
llshii et ai] (Il999l) . and iPetrovic et alJ (l2006l) ). and that a 239Mo 
helium star model would be so enormously inflated that it would 
appear as a red supergiant, if a stable envelope structure existed 
at all. This demonstrates that while the envelope physics used 
in our models leads to pushing the outer layers of stars at the 
Eddington limit to large radii, this seems to occur to a much 
lesser degree in the models of Yusof et al., with the consequence 
of a smaller mass-loss rate. 

When looking at the Eddington factor for electron scatter¬ 
ing Ee, the models of Yusof et al. behave as our models do, i.e. 
they remain mostly lower than Eg = 0.7 (cf. Sect. 13.61) . How¬ 
ever, their non-rotating SOOM© sequence is an exception, where 
Eg = 0.82 is reached at core hydro gen exhaus t ion. In comparison 
to the stellar wind calculations of Vink et aki (1201 Ih . who found 
a strong mass-loss enhancement for stars above Eg = 0.7, we 
would again expect that the mass at core hydrogen exhaustion 
of this sequence could be considerably lower than that found by 
Yusof et al. 

Our models show the largest mass-loss rates during their 
Wolf-Rayet phase (cf. Sect. 3.2). Nevertheless, it is not the 
choice of the Wolf-Rayet mass-loss prescription which makes 
our non-rotating models lose more mass than the comparable 


models of Yusof et al. We And, for example, that for our slowly 
rotating 500M© model s at the time of the max imum mass-loss 
rate, the prescription bv iNugis & Earned (120001) leads to a mass- 
loss rate for our stars which is two times larger than the one we 
use. We conclude that the high pre-Wolf-Rayet mass loss im¬ 
posed by the inflated envelopes of our models plays a crucial 
role in explaining the smaller TAMS-mass of our very massive 
slow rotators compared to those of Yusof et al. 

4.3. Rotational mixing of helium 

Another difference between the models of Yusof et al. and our 
models concerns the rotational mixing of helium. As shown in 
Sect. 13.31 the mixing of helium in our models is practically ab¬ 
sent below a threshold rotational velocity, and nearly complete 
above the threshold. The mixed models suffer more mass loss 
since they become Wolf-Rayet stars earlier. As the threshold 
velocity decreases for higher masses, this results in a bimodal 
behaviour of the total mass lost during core hydrogen burn¬ 
ing (Fig. I2TI1 . Our models with an initial rotational velocity of 
~ 300km/s follow the non-rotating models below ~ 150M©, 
while they behave in a similar way to our fastest rotators for 
higher mass. In the models of Yusof et al., the mixing of helium 
is gradually increased for larger rotational velocities. Therefore, 
their rotating models end core hydrogen burning with masses 
that exceed those of our homogeneously evolving models. 

A more detailed comparison with the models of Yusof et al. 
is difficult because they considered only two rotational veloc¬ 
ities, i.e. tirot.i = Okm/s and Vrot.i - 400km/s. This may also 
prevent a clear disc rimination of their mode l s from ours through 
observations, since iRamfrez-Agudelo et al.l (1201 3l) showed that 
only very few O stars are expected with Wrot.i = 0 km/s and Vrot.i - 
400 km/s, while the majority rotate with Wiot,i = 100... 200 km/s. 

5. Summary 

We present a detailed grid of stellar evolution models for single 
stars with initial masses from 70 to 500 M© and rotation rates 
up to 550 km/s. W e used the same physics and assumptions as 
iBrott et ^ (I2OI lah did for stars in the mass range 5-60M©, and 
this new grid therefore is an extension of their work. The ini¬ 
tial composition of our models corresponds to abundance mea¬ 
surements for massive stars in the Large Magellanic Cloud. We 
follow the evolution of the stellar models through their core hy¬ 
drogen burning phase, with some of them computed well beyond 
this stage. 

Given the hig h fraction of close binaries in massive stars 
(ISana et al.l l2012l) . we cannot hope to obtain a complete pic¬ 
ture o f massive star evolution from our models dde Mink et al.l 
1 20141), which is true in particular for the most massive stars 
(ISchneider et al.ll2014l) . On the other hand, our rotating mod¬ 
els may constitute a fair approximation for those mass gain¬ 
ers and mergers in close binaries which rejuvenate after the 
binary interaction (iBraun & Langed Il995l) . even though their 
age, their d etailed surface abund ances, and their spin may still 
be peculiar (Ide Mink et al.ll2()T^ . In addition, since most mass 
donors in post-interaction c lose binaries either remain unob¬ 
served (|de_N^k_e^l.l 1^ 1 4l) or show e x treme surface enrich¬ 
ments (lLangeill2012h . and as ISana et aki (1201 2l l2013al) suggest 
that 30-50% of the massive stars in 30 Dor appear not to have 
a close companion, a comparison of our single-star models with 
observed very massive stars will still be meaningful. 

We find that stellar rotation can influence the evolution of 
our models significantly, mostly through mixing of helium. We 


Article number, page 15 ofl21l 


















































A&A proofs: manuscript no. p9 


find two threshold initial rotational velocities, which both de¬ 
crease with increasing initial mass, one below which the rota¬ 
tional mixing of helium is negligible, and a second one above 
which the models undergo quasi-chemically homogeneous evo¬ 
lution. For initial rotational velocities in between the threshold 
values, the models start chemically homogeneously and develop 
helium-enriched surfaces, but transit to normal evolution as a re¬ 
sult of spin-down. 

Above an initial mass of ~ I6OM0, we find that quasi- 
chemically homogeneous evolution can also be achieved through 
mass loss, which in fact ensures that all our more massive models 
end core-hydrogen burning as nearly pure helium stars, indepen¬ 
dent of their initial rotation rate. We find that single stars with 
initial rotational velocities below 300 km/s need to be more mas¬ 
sive than ~ IOOM0 in order to achieve a surface helium mass 
fraction above 35% (cf. Fig.fTOli. 

Significant mixing of trace elements, including nitrogen and 
boron, is found well below the threshold velocities mentioned 
above. Because of the high core temperatures in our very mas¬ 
sive models, the NeNa- and MgAl-cycles are also activated and 
lead to surface abundance changes of the involved isotopes. A 
tripling of the surface nitrogen abundance occurs in our models 
above an initial rotational velocity of 150 km/s and a mass of 
~ 25 Mq, while a doubling at this rotation rate occurs even down 
to 6M0 (Fig.[TTJ. Because of the strong mass loss and large con¬ 
vective core masses, the surfaces of our models above ~ IOOM0 
become nitrogen-rich irrespective of their rotation rate. 

Our zero-age models above ~ 150Mo, as well as more 
evolved models down to ~ 4QMq show significant envelope in¬ 
flation because the true Eddington factor inside the stellar en¬ 
velope approaches or even exceeds the value of unity, while 
the electron-scattering Eddington factor remains below or near 
a value of Eg = 0.7 (cf. Sect. I3.6l l. Consequently, our zero-age 
main sequence models show a maximum surface temperature of 
56 000 K at about I8OM0, with lower values for higher masses 
(Fig. 13. 

During core hydrogen burning, the envelope inflation drives 
our slow rotators in the mass range 50 - IOOM0 to effective tem¬ 
peratures near or below lOOOOK. At these cool temperatures, 
partial hydrogen recombination leads to high opacities, which 
in turn produces true Eddington factors well above one in the 
envelopes of these models. While a detailed study of this phe¬ 
nomenon remains to be done, we speculate that this feature is 
related to the outbursts of luminous blue variables. 

The inflation and the corresponding decrease of the surface 
te mperature lead to e nhanced mass loss compared to the models 
of lYusof et alJ(l2013h . The enhanced mass loss has strong conse¬ 
quences for the final fate of very massive stars. In this context it 
is important to em phasize the metallicity dependence of inflation 
(llshii et alJll999h . The reduced masses (Eig. |2TJ exclude very 
massive stars of the considered metallicity as proge nitors of pair- 
instab ility supernovae. This is in agreement with iLanger et al.l 
(l2007h . who demonstrated that models with the same envelope 
physics as used here may produce pair insta bility supernoyae at a 
metallicity of Z0/2O (iKozyreya et alJ20'l4l) . This is in contrast to 
the finding of Yusof et al., who suggest that slow rotators aboye 
~ 3OOM0 end up with CO-cores aboye 6QMq at a metallicity of 
that of the EMC. 

A further consequence of the strong mass loss of our inflated 
models is their distinct spin-down during core hydrogen burn¬ 
ing (Eig. |6l). This preyents all our models aboye ~ 6QMq from 
producing a long-duration gamma-ray burst upon coll apse. Erom 
the mod els considered here and including those from iBrott et al.l 
(1201 lah . the highest chance to produce long-duration gamma-ray 


burst at EMC m etallicity from the chemically homogeneous eyo- 
lution scenario (lYoon & Langeiil2005l : IWoosleY_&jfe geill2006f) 
may occur in the mass range 20 - 3QMq (iLangei 12012 ). 

The VLT-ELAMES Tarantula Suryey (lEyans et akl 1201 ih 
proyides one of the major motiyations for this study. Our models 
will be extensiyely compared with these and other obseryations 
in the near future. To this end, we proyide the main outputs of 
our stellar eyolution models a s electronic tables ( stored at CDS) 
in the same format as those of iBrott et al.l (1201 lal) . Additionally, 
we proyide sets of isochrones in AppendixlAl 
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Appendix A: Isochrones 

The evolution of our rapidly rotating stellar models is influenced 
significantly by rotation (see Sect. 13.11) . Therefore, isochrones 
of rotating models calculated for the same age differ from 
isochrones of non-rotating models. 



Fig.A.l. Isochrones of non-rotating stellar evolution models in the mass 
range 10 - 500 Mq are depicted for ages up to 5.4 Myr. The different line 
colours explained in the figure key indicate the age of every isochrone. 
Models in lOiWo-steps are highlighted by filled circles. 

Figure lA.ll depicts isochrones of non-rotating models for 
ages up to 5.4 Myr in the HR diagram. Eight different isochrones 
are shown with age steps of 0.8 Myr. The core hydrogen burning 
stellar evolution models presented in Sect. |3] are used to gener¬ 
ate the isochrones. Because more massive stellar models have 
shorter lifetimes, older isochrones terminate at the less massive 
model at the upper end of the track. 

For a given initial composition and age, isochrones span an 
area in the HR diagram when different initi al surface rotationa l 
velocities are considered simultaneously (iBrott et al.l 1201 lal) . 
Isochrones of 16 different ages from 0.2 to 6.2 Myr for rotat¬ 
ing stellar models are shown in Fig. IA.2I and IA.3I for differ¬ 
ent ages. Switching from inhomogeneous to chemically homo¬ 
geneous evolution, the isochrone is located at higher effective 
temperatures and luminosities than for the non-chemically ho¬ 
mogeneous evolution. 

The more massive a stellar model is, the earlier it reaches the 
point in the HR diagram where the stellar evolution track starts to 
evolve blueward. The related turn in the isochrones in Fig. lA.2l is 
first visible for 1.4 Myr (bottom left panel). The minimum effec¬ 
tive temperature depends on the surface rotational velocity at the 
ZAMS, which is reflected by the isochrones of different rotation 
rates. The 2.2 Myr and 2.6 Myr isochron es show that the more 
massive a stellar model is, the earlier the iHamann et al.l (Il995h 
mass-loss rate is applied to the stellar evolution calculation. The 
isochrones therefore show a decrease in luminosity for the most 
massive models. 

Figure IA.4I depicts a population synthesis of stars with the 
age and distributions of mass and surface rotational velocity 
given in Ta ble lA. 1 1 The calcu lation was done using the code 
Starmaker (iBrott et al.l 1201 Ibl) . using the parameters listed in 
Table Al. 

It can be seen that the randomly drawn stellar models with 
the same age do not lie on one line, but instead spread over a 
certain area in the HR diagram. The initial mass distribution de¬ 
termines the model density along the line corresponding to the 
isochrone of non-rotating stellar models of this age. We choose 


Table A.l. Parameters used in our population synthesis calculation (cf. 
Fig.lAAl. 


parameter 

age 

1.5 Myr 

velocity distribution 

Gaussian distribution IBrott et al. 201 Ibf 


(cr = 141 km/s,/j=100km/s) 

velocity range 

0 - 500km/s 

mass distribution 

uniform distribution 

mass range 

10-500 Mo 


a flat mass distribution for the simulation to have a better view 
of the behaviour of the most massive stars. The velocity distri¬ 
bution on the other hand shapes the stellar model density as a 
function of the effective temperature for (roughly) constant lu¬ 
minosity. The colour coding indicates the number of stars within 
one pixel of 500 K and log(L/Lo)=0.05. 



Fig. A.4. STARMAKER population-synthesis calculation for 1.5 Myr (and 
other parameters as explained in Table lA. fl . The luminosity is depicted 
as a function of the effective temperature. The colour coding explained 
at the right bar indicates the number of stars within one pixel of 500 K 
and log(L/Lo)=0.05. 

Figure IA.4I shows similar information as discussed in 
Fig. lA.2l Additionally, it gives the probability of observing a star 
for the conditions given in Table lA.ll depending on the surface 
rotational velocity distribution. 

Isochrones for slow rotating stellar evolution models that 
undergo non-chemically homogeneous evolution do not dif¬ 
fer significantly in the HR diagram. The probability of ob¬ 
serving a star along these isochrones is highest. The change 
from non-chemically homogeneous to quasi-chemically homo¬ 
geneous evolution occurs at a small velocity range related to a 
significant change of the position in the HR diagram.Therefore, 
the probability of observing a star between isochrones of slow- 
rotating and rapidly rotating models is small. Again, isochrones 
of rapidly rotating stellar evolution models are close in the HR 
diagram, leading to a higher probability of observation. 
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Fig. A.2. Several isochrones of different ages and rotating stellar evolution models are shown in the HR diagram. The surface rotational velocity 
at the ZAMS is chosen in steps of lOOkm/s from non-rotating to 400km/s and additionally 450km/s. Three initial surface rotational velocities are 
highlighted in particular. The isochrones of non-rotating models are shown in red, 400 km/s in blue, and 450 km/s in orange. 
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Fig. A.3. Several isochrones of different ages and rotating stellar evolution models are shown in the HR diagram. The surface rotational velocity 
at the ZAMS is chosen in steps of lOOkm/s from non-rotating to 400km/s and additionally 450km/s. Three initial surface rotational velocities are 
highlighted. The isochrones of non-rotating models are shown in red, 400 km/s in blue, and 450 km/s in orange. 
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Table B.l. Summary of important stellar characteristics. 






